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Abstract 


We present in this paper the current status of searches for leptoquarks at 
D0. These results include the use of next leading order theoretical predictions 
for the cross section for pair production of leptoquarks at hadron colliders. 
We also present a new optimized analysis for first generation leptoquarks 
with significant increase in sensitivity relative to earlier searches using D0 
data. The mass limits derived from this first generation leptoquark search 
are relevant to the recently reported high mass events at HERA. 


I. INTRODUCTION 

Leptoquarks are hypothetical exotic particles with both lepton and color quantum num¬ 
bers; they are color triplet bosons with fractional charge. They are produced in pairs at 
hadron colliders via strong interactions. A leptoquark will decay via an unknown coupling, 
A, to a lepton and a quark. The production of leptoquark pairs in pp collisions is insensitive 
to A, and we are not concerned with this coupling as long as it is greater than 10“^^; other¬ 
wise, the leptoquarks will not decay in our detector and escape detection. We assume that 
leptoquarks occur in generations; this is to say that leptoquarks of one generation couple 
exclusively to leptons and quarks of the same generation. For example, a hrst generation 
leptoquark will only couple to electrons, electron neutrinos and u or d quarks. With this 
assumption of generations, we are able to bypass strict bounds (> 1 TeV) on leptoquark 
masses derived from limits on flavor changing neutral currents. 

Since leptoquarks are produced dominantly in pairs at hadron colliders, their signature 
would be two leptons plus two or more jets. A free parameter of the model for leptoquarks 
is the branching fraction of the leptoquark to charged lepton plus quark, /?. For f3 = 1.0, 
leptoquarks decay 100% to a charged lepton plus quark. In the case for a hrst generation 
leptoquark pair, the signature would be two isolated electrons (e'''e“) plus at least two 
jets. Additional jets could arise from initial state radiation or hnal state radiation. Since 
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we assume leptoquarks do not couple outside their own generation, for first generation 
leptoquarks, we would not expect to see muons or evidence of tans in the signature. 

Some published and previous preliminary conference limits on leptoquarks masses are 
shown in tables | and |1[ These give a brief history of mass limits from the Tevatron. The 
previously published limits from CDF and D0 have used Leading Order (LO) theory. Next 
Leading Order (NLO) theory predicts cross sections that are 30-50% higher than LO theory. 
NLO was used for the CDF limits at recent conferences []^ [Q and will be used in the newest 
limits presented in this paper from D0. 


II. FIRST GENERATION LEPTOQUARKS 


DO searches for hrst generation leptoquarks primarily in the dielectron plus two jet and 
electron plus missing transverse energy ) plus two jet signatures corresponding to the 
decay modes eeqq and evqq. The search for hrst generation leptoquarks in the dieletron plus 
two jet signature has been greatly modihed since the 1997 Moriond Conferences 0. After 
the announcement of the high mass events from HERA [10|, DO has retooled the analysis 


to increase the acceptance for leptoquarks by increasing the acceptance for electrons and 
by optimizing the event section for higher mass leptoquarks. The data sample used in this 
analysis is the full Run I data sample of 123 pb“L The single electron plus plus two jets 
signature search is also optimized at higher masses. We are also using the newly available 
NLO theory predictions for the production cross section. 

The electron identihcation for one of the two electrons in the dielectron search has been 
loosened to increase the acceptance for leptoquarks. For the loose selection, the track re¬ 
quirement is removed. This has effectively doubled the acceptance for dielectrons. The 
relaxing of the track requirement for the one electron does not signihcantly increase the 
QCD background relative to the tight-tight electron selection. 

The basic eejj event selection is given in the following list: 

• two electrons with Et > 20 GeV with | r; |< 1.1 or 1.5 <| r/ |< 2.5, 

• two jets with Et > 15 GeV with | ?; |< 2.5 where the jets are reconstructed with a 

cone algorithm using a cone of radius R = y/Ar]‘^ -|- = 0.7, 

• exclude events with dielectron invariant masses between 82 to 100 GeV/c^. 

The dielectron mass cut is a three sigma cut. With this basic selection, we have 101 events 
in our candidate sample. The main backgrounds come from Drell-Yan, QGD with jets faking 
electrons, and Top. We estimate 66.8±13.4 events for Drell-Yan, 24.3±3.6 events from QGD, 
and 1.8±0.7 events from Top for a total of 92.8±13.8 background events. 

Further event selection is optimized for leptoquark masses of 200 GeV/c^ and greater. 
The optimization is based on signal and background Monte Garlo event samples and on 
QGD data event samples. Several variables have been studied for the optimization; variables 
involving sums of energies and transverse energies {Et) of electrons and jets, variables of 
event shape, variables of reconstructed invariant masses, and variables involving constraints 
with leptoquark masses have been used. Systematic grid search |]^ and neural network 
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techniques are used for the optimization studies. We approached the optimization in two 
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ways. The first is a discovery search where we optimized signal (S') over the square root 
of the background (B), S/^/B. The second is for setting a limit where we maximized the 
signal for a given target background of about 0.4 events. This level of background gives a 
probability of about 70% of seeing 0 background events. 

We hnd that a transverse energy variable dehned as St = + E^ is the most 

jets 

optimal variable. Here, the jets are required to have Et > 15 GeV to be included in the sum. 
Figure |I| shows the output of an example grid search. Here is plotted the number of signal 
events verses the number of background events. Each point in the plot represents a different 
set of cuts. The electron and jet E^s are fixed at the basic cut level as given in the hgure. 
The upper line of dots are for the St cut. The lower line of dots is for a leptoquark mass 
constrained variable dehned as: DM/M (200) = ^{Meji — M'^qY + {Mej 2 — )2/M|q . 

Mgji and Mf.j 2 are the invariant masses of the two combinations of the electrons and jets 
such that the difference in the two masses is minimized, and is the mass of a 200 

GeV/c^ mass leptoquark in this example. The wide band of dots is the combination of both 
variables used in the selection. From this plot we see that that the St variable does better 
than the mass constraint variable at our target background of 0.4 events. We hnd that a 
St > 350 GeV gives a background of about 0.4 events. None of the 101 candidate events 
pass this cut. The diherential and integrated St distributions for the data and background 
predictions are given in Figs. and ^ The background distribution matches well the data 
distribution. 

We have also studied the mass properties of the 101 dielectron candidate event sample. 
We used a 3G (3 constraint) ht where we require balance of the transverse energy in the 
event and we require the masses of the two reconstructed electron-jet systems be equal. The 
combination of electron with jet is chosen such that the diherence in masses is minimized 
before the htting. The 3G mass ht for a 225 GeV/c^ mass leptoquark sample {j3 = 1.0) can 
be seen in Fig The peak of the distribution is about 10% low but the mass resolution is 
good at about 15 GeV/c^. In Fig. ^ is the 3G htted mass for the data shown as the points 
with errors. The background prediction is shown as the solid line histogram. The data and 
background prediction agree well. The dashed line histogram is the 3G htted mass for a 200 
GeV/c^ leptoquark sample. We expect about 6 events from a 200 GeV/c^ mass leptoquark. 
By comparison we expect about 20 events from a 160 GeV/c^ mass leptoquark and 3 events 
from a 225 GeV/c^ mass leptoquark. All histograms are normalized to 123 pb“^. 

One might note the two events at the high end of the mass distribution in Fig. These 
two events have values of St that are signihcantly lower than our cut of 350 GeV. This can 
be seen in Fig. Here we plot the St variable versus the 3G htted mass for the background 
prediction (upper left), for a 225 GeV/c^ mass leptoquark sample (upper right), and for 
the data (lower left). The two events in question have S'^’s less than about 200 GeV. From 
Fig. ^ we see that the data and background prediction are again very similar. 

Given that we see no dielectron leptoquark signal in our data, we can proceed to set a 
limit on the production cross section for leptoquark pairs, and by comparing to theory we 
can set limits on the leptoquark mass. The efficiency for detecting the dielectron plus two 
or more jet signature from leptoquark pair production as a function of leptoquark mass is 
given in Fig. |^. The lowest curve is the total efficiency. The errors bars represent the total 
statistical plus systematic uncertainties of about 13%. These uncertainties are listed here: 
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• energy scale: 2-5% 

• electron identification: 5% 

• acceptance: 5% 

• gluon radiation: 7% 

• parton distribution functions and Q^: 7% 

• luminosity: 5% 

• Monte Carlo statistics: 2% 


Given the efficiency and uncertainties, we calculate |]^ a 95% CL limit on the production 
cross section times as a function of leptoquark mass. This is given in Fig. Also given 
is the NLO theoretical prediction. The band represents the range in cross section prediction 
as the renormalization scale changes from one half the leptoquark mass (upper boundary: 
/i^ = 1/4 X MIq) to twice the leptoquark mass (lower boundary: = 4 x MIq). The 

intersection of our experimental limit on the cross section and the lower boundary of the 
theory prediction gives a 95% CL limit of 225 GeV/c^ on the mass of the leptoquark for 


(3 = 1.0. 

We have also searched for first generation leptoquarks with the single electron plus 
plus two or more jets signature. Recall that for this signature one leptoquark decays to a 
electron plus quark, and the other leptoquark decays to electron neutrino plus quark. In 
this case we are most sensitive to the branching fraction (3 = 0.5. We have 103.7 pb“^ for 
this search. The basic event selection is given as follows: 


• one electron with Ex > 25 GeV, \ rj \< 1.2, 

• two jets with Ex > 25 GeV, \ t] \< 1.0; a third jet with Ex > 25 GeV is allowed 
in the event; additional jets with Ex < 25 GeV are not used to veto events; jets are 
reconstructed with cone algorithm (R = 0.7), 

• $T > 40 GeV, jets and $x must have A(/) > 0.25 and | tt — A0 |> 0.25, 

• Sx > 170 GeV where Sx = Y^E^^ + Ex'-, jets must have Ex > than 15 GeV to be 
included in sum, 

• exclude events with isolated muons, since muons are not part of the first generation 
leptoquark signature. 


With the basic selection we find 32 candidate events. The background sources to this 
signature are W boson -|- 2 jets, top, and QCD. The background estimates from these 
sources are 19.6±4.2 events from W boson -|- 2 jets, 9.0±2.7 events from top, and 1.1±0.4 
from QGD for a total of 29.8±5.0 events. 

The transverse mass of the electron and li^x is given in Figs. ^ and 


10. The data 


is given as the points with error bars representing the statistical uncertainties, and the 
background prediction is the dashed histogram. We see very good agreement between the 
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data and background predictions before and after the Ht cut. The arrow in the hgures 
represents a cut of 100 GeV on . This cut eliminates all but one event. This event has 
an electron plus four jet topology and is in fact a top candidate. This is supported by the 
background estimates of 0.52±0.28 events from W boson + 2jets, 1.55±0.48 events from 
top, and 0.41±0.41±0.20 events from QCD for a total background prediction of 2.5±0.6 
events. 

The last component of the event selection is based on a variable, 5M dehned as 5M =\ 
Mej — Mlqi I /Mlqi. Mlqi is the mass of the hrst generation leptoquark, and is the 
electron - jet invariant mass. The jet that is chosen from the two or more jets in the event 
minimizes 6M. Optimizations of this variable using signal and background Monte Carlos 
reveal a cut of SM < 0.2 is optimal. We hnd that there are no events left in our sample for 
leptoquark masses greater than 140 GeV/c^. 

With the signal efficiencies shown in Fig. we calculate the 95% CL limit on the 
production cross section times 2/3(1 — /3). This is given in Fig. as function of leptoquark 
mass. The theoretical prediction is the NLO production cross section times 2/3(1 — /3) for 
[3 = 0.5. Our 95% CL limit drops at 140 GeV/c^ because for masses greater than 140 
GeV/c^ our one remaining event does not survive the 5M cut. The intersection of our limit 
with the lower theoretical prediction = 4 x M|q) gives a mass limit on hrst generation 
leptoquarks of 158 GeV/c^ for (3 = 0.5. 


III. SECOND GENERATION LEPTOQUARKS 


The search for second generation leptoquarks involves signatures of muons and jets. The 
search described here is for two isolated muons plus two or more jets. This search uses 94.4 
pb“^ of data. The basic event selection is the same as the standard D0 dimuon top quark 
selection [|T^ which because no direct cut is used also does well for the leptoquark dimuon 
event selection. The major backgrounds to second generation leptoquarks are Drell-Yan and 
top. The event selection is listed here: 


• two isolated {R{fi,jet) > 0.5) muons with Et > 15 GeV and \ rj \< 1.0, 

• two jets with Et > 20 GeV and | ?7 |< 2.5, 

• < 160° if I + r/^2 |< 0.5, 

• the dimuon invariant mass, greater than 10 GeV, 

• HT{jets) > 100 GeV, 

• Z-kinematic ht probability, p(x^), less than 1%. 


Here Ht is the sum of the Et of jets for jets with Et > 15 GeV. The back-to-back cut in 0 
on the muons is for cosmic muon rejection. The Ht cut is intended to signihcantly reduce 
the Drell-Yan background. This can be seen in Fig. ^ where we show the Ht distribution 
for Drell-Yan, top, and leptoquark (M = 160 GeV/c^) Monte Garlos. We see that the Ht 
cut of 100 GeV rejects a large portion of the Drell-Yan [Z /i/i) background. It keeps most 
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of the top and leptoquark events. With these cuts we have one event with an estimated 
background of 0.97±0.20 events from top and Drell-Yan. 

To reduce the top background we consider the 0 distribution of the two muons and the 
two highest Et jets in the candidate events. Top pair (tt) events producing the dimuon 
(isolated muons) signature contain $rp , the source of which is two neutrinos from the decay 
of the two W bosons that came from the decay of the tt pair. The dimuon signature from 
second generation leptoquark pair production has no from real sources like neutrinos. 
It is possible that the neutrinos in the top events will concentrate in one 0 region of the 
detector causing the jets and muons to concentrate on the opposite side to balance the 
event. This will tend to produce a large gap in 0 between two of the muons or jets. One 
can imagine a pie cut up into four pieces by the two muon and two jet directions. Using 
this analogy, a top event could cut out a piece of the pie that is more than half the total 
pie. Leptoquark events would tend to cut the pie into more equal pieces. The distribution 
of the maximum 0 gap is given in Fig. Ti for leptoquark, top, and Drell-Yan Monte Carlo 
samples. If we require that the maximum 0 gap be less than 180°, we reject a signihcant 
portion of the top background while we retain nearly all of the leptoquark signal. 

With the 0-gap cut we have no events left in our sample. The last event is actually a tt 
candidate with a very large 0-gap. The total efficiencies for the signal varied from 3.8%±0.5% 
to 12.6%±1.3% for second generation leptoquark mass that varied from 100 GeV/c^ to 260 
GeV/c^. The total statistical and systematic errors are 10-15%. We calculate the 95% CL 
limit on the production cross section time 0^ for second generation leptoquark pairs and 
compare this to the NLO theory = M|q). In Fig. ^ we show the preliminary 95% CL 
limit exclusion contour (single line) of 0 vs second generation leptoquark mass. For 0 = 1.0 
we have a mass limit of 184 GeV/c^, and for 0 = 0.5 we have a limit of 140 GeV/c^. The 
hatched regions are the limits from a previous analysis using only the 1992-1993 (lA) 
data set. 


IV. THIRD GENERATION LEPTOQUARKS 

We’ve searched for signatures of third generation leptoquarks of charge 1/3. These 
charge 1/3 leptoquarks decay to b quark plus z/,-. Our search reach for third generation 
scalar leptoquarks is below the top mass of about 170 GeV/c^, so we assume that third 
generation scalar leptoquarks of charge 1/3 decay 100% of the time to b quark plus 

In the event selection we require $rp > 35 GeV, two jets at least one of which has a muon 
tag, and topological cuts. The untagged jets are required to have Et > 25 GeV, and the 
tagged jets are required to have Et > 10 GeV (excluding the muon Et). We have total 
efficiencies of 2-5% for third generation leptoquark masses between 100 GeV/c^ and 300 
GeV/c^. The major background sources to this signature are top, W and Z bosons plus two 
jets, and QGD multijets. For these cuts we see two events in the full Run 1 data sample 
(about 20 pb“^ for our selected trigger) with an expected total background of 3.1Y0.9 events. 

In Fig. ^ we show the 95% GL limit on the cross section times (1 — 0)^ as the stars 
connected by the dotted line as a function of third generation leptoquark mass. The solid 
line is the NLO theory (p^ = M^q) for scalar leptoquarks and the dashed line is the LO 
theory for vector leptoquarks with Yang-Mills coupling. From this plot we see that we set 
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a limit on the mass of scalar third generation leptoquarks (Q = 1/3) of 98 GeV/c^, and we 
set a limit on third generation leptoquarks of 201 GeV/c^ (/3 = 0). 


V. CONCLUSION 


We have searched for three generations of leptoquarks with diagonal couplings to leptons 
and quarks. We have found no evidence of a leptoquark signature in the D0 Run I data. A 
summary of the preliminary mass limits for the three generation of leptoquarks is given in 


table III 


For the future, we expect that our mass reach for lower (3 in the second generation 
search to improve when we add the search for the single muon plus plus jets signature 
to this analysis. We expect to greatly improve our search reach in both the hrst and second 
generation leptoquark searches at low [3 when we have added our $rp plus jets searches to 
these analyses. Finally, we will also have limits for vector leptoquarks. 
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TABLE 1. 

colliders. 

TABLES 

Brief history of published 95% CL mass limits 

on scalar leptoquarks from hadron 

Experiment 

Signature 

fj 

95% CL mass limit (GeV/c^) 

CDF i 

eejj 

1.0 

113 

CDF i 

eejj 

0.5 

80 

D0 § 

eejj, eujj 

1.0 

130 

D0 i 

eejj, eujj 

0.5 

116 

HERA (HI) g 

ej 

1.0 

2750 

CDF i 

mj 

1.0 

131 

CDF § 

mj 

0.5 

96 

D0 § 

f^^j, ^^^jj 

1.0 

119 

D0 i 


0.5 

97 

CDF i 

rrjj 

0=4/3,2/30 

99 


TABLE H. Brief history of 95% CL 
recently presented at conferences. 

mass limits on 

scalar leptoquarks from hadron colliders 

Experiment 

Signature 

P 

95% CL mass limit (GeV/c^) 

D0 0 

eejj, eujj 

1.0 

175 

D0 0 

eejj, eujj 

0.5 

147 

D0 0 

eejj, eujj 

0.0 

81 

CDF i 

eejj 

1.0 

210 (NLO)0 

D0 0 

mj, IJ-i'jj 

1.0 

167 

CDF 0 

mj 

1.0 

197 (NLO) 

D0 0 

bbu-rUr 

Q = 1/3 

80 

CDF i 

TTjj 

Q=4/3,2/3 

no (NLO) 


^This limit is sensitive to the leptoquark - lepton - quark coupling; limits from HERA assume 
that A = o; 0 ni‘ 

^ For this third generation leptoquark limit the Top quarks is not relevant, so (3 is not relevant. 

^This limit from CDF uses Next Leading Order (NLO) theory to set limits on the leptoquark 
mass. 
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TABLE III. 

Summary of scalar leptoquark 

mass limits. 

Generation 

P 

95% CL limit (GeV/c2) 

comment 

Irst 

1.0 

225 

NLO (^2 = ^mIq) 

Irst 

0.5 

195 

NLO (^2 = 4m2q) 

2nd 

1.0 

184 

NLO (^2 = mIq) 

2nd 

0.5 

140 

NLO (;u2 = mIq) 

3rd 

- 

98 

Q = 1/3, NLO (/i2 = mIq) 
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FIG. 1. Output of grid search. We can see by this that the St variable alone optimizes the signal 
to background ratio better than the mass constrained variable, DM/M(200), and the combination 
of both. 
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FIG. 2. The differential distribution of St for data and background predictions. The basic 
event selection is applied. 
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FIG. 3. The integrated distribution of St for data and background predictions. The basic event 
selection is applied. 
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FIG. 4. The 3C mass fit for a 225 GeV/c^ mass leptoquark sample, f3 = 1.0. 



FIG. 5. The 3G mass fit for the data (points) with the basic event selection shown with back¬ 
ground prediction (solid line historgram) and 200 GeV/c^ mass leptoquark sample (dashed line 
histogram). Preliminary. 
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FIG. 6. The St vs 3C mass fit for the data with the basic event selection, background prediction, 
and 225 GeV/c^ mass leptoquark sample. 
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FIG. 7. The efficiency for detection of the dielectron plus two or more jets leptoquark signature. 



FIG. 8. The 95% GL limit on the production cross section times /9^. Also shown as the band 
is the NLO theoretical prediction Q. 
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M^(ev) GEV 

FIG. 9. The e — v transverse mass distribution prior to the Ht cut. The number of data events 
is 69. Data are the crosses, and the dashed histogram is the background prediction. The arrow 
indicates the mass cut. 



M.j,(ev) GEV 

FIG. 10. The e — v transverse mass distribution after to the Ht cut. The number of data 
events is 32. Data are the crosses, and the dashed histogram is the background prediction. The 
arrow indicates the mass cut. 
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FIG. 11. The efficiency for detection of the electron plus plus two or more jets leptoquark 
signature. 



FIG. 12. The 95% CL limit on the production cross section times 2/3(1 — /3) for first generation 
leptoquarks. Also shown is the NLO theoretical prediction. 
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FIG. 13. The Ht distribution for top and Drell-Yan backgrounds and a 160 GeV/c^ mass 
leptoquark signal sample. Our cut is indicated by the arrow. 
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FIG. 14. The <^-gap distribution for top and Drell-Yan backgrounds and a 160 GeV/c^ mass 
leptoquark signal sample. Our cut is indicated by the arrow. 


18 
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FIG. 15. The 95% CL limit exclusion contour for second generation leptoquarks. Plotted is the 


branching fraction, (5, vs leptoquark mass. We exclude the region to the left of the curve. We use 
NLO theory = M£q) to determine this contour. 


3rd Generation, Charge 1 /3 LQ Pair Production from DO 
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FIG. 16. The 95% CL limit on the production cross section for charge 1/3 third generation 
leptoquark. We show here as the solid line the NLO theory = M/q) for the scalar leptoquarks. 
The dashed line is LO theory for vector leptoquarks. 
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